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Abstract 

A simple counterexample against the Vlasov equation is put for- 
ward, in which a magnetized plasma is perturbed by an electromag- 
netic standing wave. 

In today's physics, the Vlasov equation serves as the very foundation for 
a great number of formulations and calculations in fields ranging from star 
physics to controlled fusion. 

The equation states that the distribution function f(t, r, v) of a collection 
of moving charged particles (called a plasma in the usual physics language) 
satisfies [H [2] 

f + v .§f + JL (E + vxB ).§Uo, (!) 

ot or m av 

where (E, B) represents the electromagnetic field experienced by the parti- 
cles of charge q and mass m at the position r and with the velocity v. It 
is sometimes called the collisionless Boltzmann equation in statistical me- 
chanics due to the fact that the equation takes almost the same form of 
the ordinary Boltzmann equation with the collision term neglected. In an 
alternative form, Eq. (pQ) is written as[TJ [3] 



df[t,r(t),v(t)} 



dt 



= 0, (2) 
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where the subscripts imply that the derivative is performed along a par- 
ticle's path. This expression can be interpreted as the path-invariance of 
the distribution function in the six-dimensional phase space spanned by the 
position r = (x,y,z) and velocity v = (v x ,v y ,v z ). Recognized as a good 
statistic approximation by this community, Eq. ([T]), or ([2D, is very often 
used to investigate the dynamical behavior of dilute plasma. 

Being exposed by the fact that the regular Boltzmann equation suffers 
substantial difficulties [H [5], the author has somehow been compelled to find 
out in what situations the Vlasov equation can be falsified. It turns out that 
a number of counterexamples against the Vlasov equation can be constructed 
in plain manners. In this paper, we shall focus ourselves on one of them, in 
which the six-dimensional physics is reduced to a one-dimensional one and 
can be understood with conceptual clarity. 

Let's start with a plasma having a strong constant magnetic field Bq 
along the z-direction. Suppose that its distribution function / = f(t, r, v) 
is initially uniform in the position space, and that it is initially independent 
of velocity when the speeds of the concerned particles are below a certain 
limit vq. Namely, we have, with v < vq and at t = 0, 

W = V = df = df_ = dl = dl = Q (3) 
dx dy dz dv x dv y dv z 

or, in another form, 

/(t,r,v)| t=0 = / c , (4) 

where f c represents a constant. For further simplicity, we also assume that 
the particles whose speeds are beyond vo occupy only a small proportion of 
the plasma and can be disregarded in our consideration. 

Adopting all the initial conditions mentioned above, we see that no mat- 
ter what kind of electromagnetic field is applied, Eq. ([1]) leads us to nothing 
but 

df 



at 



= 0, (5) 

=0 

which means that the distribution function of the system does not change 
at the beginning. More than that, Eq. ([2]) informs us that if an observer 
adheres to any moving particle, he/she finds that / is equal to f c constantly, 
irrespective of the position and velocity, which literally infers that the dis- 
tribution function remains unchanged at the beginning and the later time; 
or, in terms of macroscopic quantities there is no density change and no 
current everywhere once for all (readers may also refer to Chapter 8 of Ref. 
1, in which the linearized Vlasov equation for uniformly magnetized plasmas 
yields the same conclusion). 
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However, physical inspection tells us a very different story: applying 
many types of electromagnetic fields will change the plasma significantly 
For this to be seen in an intuitive and heuristic manner, let the magnetized 
plasma be perturbed by an electromagnetic standing wave in the form: 

Ei v = — cos(u)t) sm(k x x), Biz = sin(atf) cos(k x x), (6) 

k x u 

where a is a small constant. Note that this perturbative wave obeys Maxwell's 
equation V x E = —dH/dt. 

As well known in physics, while all particles in the plasma make gyration 
due to the existence of Bq, the gyro-center of every particle (guiding center) 
involves two types of relatively slow drifts due to the existence of Ei y and 
Bi z . The first slow drift is the E x B drift along the x-direction: 

Vx = ^ET = T~~~ET cos(wt) sin(k x x), (7) 

-DO %-DO 

and the second slow drift the grad-B drift along the y-direction 
v±p dB lz ak x v ± p 



2B dx 2ujB 



sin(o;t) sin(k x x). (8) 



Equations (JT]) and ([8]) hold when uj is much smaller than the gyro-frequency 
and the perturbed fields are relatively weak. 

If observing the plasma in terms of measuring velocity and speed, two 
things can be found. The first is that expression ([7]) apparently represents 
an average velocity along the x-direction (forming a macroscopic current if 
the plasma is a single-component one). The second is that there is a speed 
change related to every particle due to the combined effect of the grad-B 
drift and the electric field E\ y although the change, of a 2 -order in magnitude, 
is relatively unimportant in a perturbative approach. 

If the density change is of concern, something against the standard theory 
will also surface. While the grad-B drift expressed by ([8|) is uniform along the 
y-direction, causing no density change anywhere, the ExB drift expressed by 
([7]) is nonuniform along the x-direction and must result in a density change 
along the x-direction. By assuming that the gyro-radii of all the concerned 
particles are relatively small (much smaller than 2ir/k x ) and noticing the 
fact that the ExB drift is velocity- independent, the density change can 
be evaluated with ease. With help of the linearized continuity equation of 
ordinary fluid, we find that 

dni dv x n a 

= ~ n o~^~ = — 5~ cos(ujt) cos(k x x), (9) 

Ot OX -Do 
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where no is the initial density of the plasma. Fig. 1 demonstrates that 
particles in this plasma, regardless of their velocities, converge and diverge 
in the pattern revealed by ([9]). 
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Fig. 1: The E x B drift of gyrating particles along the re- 
direction caused by the electromagnetic standing wave. 

It is at this stage intriguing and essential to look at what happens in 
the phase space spanned by x, y, z, v x , v y , v z . It should be pointed out that 
this six-dimensional space is a truly abstract and weird one, serving as a 
great source of all confusions and mistakes. Investigating the change of 
/ in the phase space challenges our conceptual and experimental wisdom 
though it is surely in the domain of classical mechanics. Referring to Fig. 1, 
let's imagine that there is a small virtual box adhering to a moving particle, 
represented by AxAyAz symbolically, and imagine that we are able to count 
the number of the inside particles whose velocities are within a small range, 
represented by Av x Av y Av z symbolically. What will be seen by us when the 
box AxAyAz drifts leftwards and rightwards in the circumstances? At the 
two ending points, the density in the box AxAyAz will, as revealed above, 
change from the lowest to the highest (or vice versa) while the velocities of 
the particles inside the box keep almost unchanged (noticing that the drift 
velocity at the ending points are zero and the speed changes of the particles 
are of second order). All this means that / in the phase volume element 
AxAyAzAv x Av y Av z is not invariant along a particle's path. 

Many theoretical and practical questions then arise. If the Vlasov equa- 
tion is indeed unsound, how come there exist so many academical proofs 
showing otherwise? If the standard theory fails to yield good predictions 
about density's changes and drift currents, why haven't present computa- 
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tional simulations revealed any discrepancy between the theory and reality? 
In what way can we correctly formulate the behavior of plasma? Extensive 
effort has been made by the author to answer these questions and to find 
out new things[5j [6], but the greatest challenge was, and still is, to convince 
the mainstream of this community to have open mind on related issues. 

The author thanks professors Oliver Penrose and Keying Guan for stim- 
ulating discussions. 
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